A simple algorithm for automatic Feynman diagram generation 

Bo Xiao 1 , Hao Wang 2 *, and Shou-hua Zhu 3,4 
1 Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China 

2 Department of Astronomy, Beijing Normal University, Beijing 100875, China 
3 Institute of Theoretical Physics & State Key Laboratory of Nuclear Physics and Technology, 
Peking University, Beijing 100871, P. R. China 
O ; 4 Center for High Energy Physics, Peking University, Beijing 100871, P. R. China 

a? 

C/3 ■ (Dated: September 6, 2012) 

in! 

^ Abstract 

9* 

An algorithm for the automatic Feynman diagram (FD) generation is presented in this 
— ■ paper. The algorithm starts directly from the definition formula of FD, and is simple in 

concept and easy for coding. The symmetry factor for each FD is naturally generated. 
It is expected to bring convenience for the researchers who are studying new calculation 
techniques or making new calculation tools and for the researchers who are studying 
effective field theory. A C-program made from the algorithm is also presented, which is 
short, fast, yet very general purpose: it receives arbitrary user defined model and arbitrary 
process as input and generates FD's at any order. 
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I. INTRODUCTION 

Cross section calculation is an important and hard task in high energy physics. The 
calculation technique in high energy physics is based on the Lorentz invariant perturbative 
expansion of the S-matrix developed by Feynman, Schwinger, and Tomonaga. In this tech- 
nique, a calculation begins with generating all the FD's. The number of Feynman diagrams 
for a normal scattering process can be hundreds or thousands. To obtain all of these dia- 
grams manually would be exhausting and prone to errors. The Calculations of these FD's 
are even more tiring. Thus, many software's have appeared which can manage these steps 
automatically (see |l| for a good review). 

The rules for FD generation are definite, hence are quite suitable for compute automation. 



There are we,, known general — FD generators such as FeynArts ft H , QGRAF Q 
and the grc part of GRACE [5]. FeynArts and QGRAF generate FD's with the steps of: 1. 
generating all the connected topologies; 2. eliminating equivalent topologies generated; 3. 
inserting physical fields into the topologies to generate all the FD's; 4. eliminating equivalent 
FD's generated. The grc part of GRACE generates FD's by: 1. uses an iteratively connecting 
method to generate all the FD's; 2. eliminate equivalent FD's generated. In all these three 
generators, the post eliminating procedure is time consuming and full of tricks. There are 
also many more not so general FD generators Q-Q. MadGraph Q and AMEGIC Q adopt 



a topo 



m 



11 



ogy-based algorithm. They are restricted to the tree diagrams. The FD generators 



13] are novel and fast, but have limited functionalities and are designed for some 
special purposes. Overall, the existing algorithms for general purpose FD generation are 
some complicated that are not easy for other researchers to follow to write their own FD 
generators. 

A concept-simple and coding-easy FD algorithm could prove to be helpful to, say, re- 
searchers developing new calculation tools. For example, some researchers have being devel- 
oping new techniques of doing the loop integration in a numerical way [yj]; using self-made 
FD generators would make the process of such studies more convenient. 

The rest of the paper is organized as follows. In section [Til the simple FD generation 
algorithm is described; In section IIII} a small C-program that realize this algorithm is de- 
scribed; In section ITV] the correctness of the algorithm/ C-program is checked; In section \V\ 
conclusions are given. 
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II. THE ALGORITHM FOR FEYNMAN DIAGRAM GENERATION 



The princi ples for FD generation is contained in the well-known perturbative formula for 



the S-matrix 



Sp^ a' 1 n' 1 ;p' 2 a! 2 n' 2 ;---piaini;p2(T2n2; 



= E nT~ I d 4 Xi---d 4 x N (%---a(p' 2 a' 2 n' 2 )a(p' 1 a[n' 1 ) (1) 
xT {H{x x ) ■ ■ ■ H(x N )} a^(p 1 a 1 n 1 )a' l (p2a2n2) ■ ■ ■ $ ) • 

It's convenient to convert (OQ) into an abbreviate form, since one concerns only the FD 
generation here 

SiV.i 2 ... ~ ^ {■ ■ ■ a'Am^ ■ ■ ■ H N }{a\at ■■■} (2) 

It is also convenient to turn all the particles in the "in" state into its antiparticles in the 
"out" state 

S 12-n ~ jyj {a-nO-n-i • ■ ■ a c 2 a{}T {H^ ■ ■ ■ H N } , , (3) 

where the primes have been discarded and subscripts have been modified. This conversion 
do not change the structure of a FD. Eq. ([3]) is the base formula for the following discussion. 

Let's consider a FD generation example for the process ei<32 — > e^e^ and the number of 
interaction vertex N = 2 and the model 



H oc V^ e n (4) 

In this example, the expression (|3]) can be put into an explicit form as shown in Eq. (jSJ) 
(note: e\ and e-2 in the "in" state have been turned into their antiparticles e\ and €2 in the 
"out" state in Eq. ©) 

5 1234 ~ 2! ^ 4e3e2 ^} T i [folfafa] } (5) 

All the required FD's are contained in expression (JSJ), that is, they correspond to all the 
combinations of grouping the annihilation operators and the field operators according to the 
pairing 

(#), (eV>), (ruf>), O/O, and (#). 
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To find all the combinations, we define an order for the pairing procedure 
Pairing order: 

• Firstly, the external particles are scanned one by one from right to left (from the ones 
with smaller subscript to the ones with larger subscript) to pair with the fields in the 
T{} brace; 

• Secondly, after all the external particles are scanned, the remaining unpaired fields in 
the T{} brace are scanned one by one from left to right (from smaller subscript to 
larger subscript) to pair with other fields 

• When an operator (an external particle or a field in the T{} brace) is pairing with the 
fields in the T{} brace, it scans the fields from left to right. 

Following this pairing procedure we got for Eq. (JSJ) four combinations, they are put in 
Eq. ([6]) in the order of their production 

±T { [(ea&XexVi)] [{e^WM^M] } 
JjT { [(es^Xe^i)] [(e 2 ^ 2 )(e 4 ^ 2 )(0 1 2 )] } 

±T { [(e 2 &Xe 4 Vi)] [(es^Xei^X^)] } 
±T { [(e 3 ^i)(e#i)] [M2)(eM{<t>M] } 

The four combinations in Eq. (jSJ) correspond to four diagrams in Figure [I] respectively. 




FIG. 1: four Feynman diagrams corresponding to the four combinations in (6). 

In Figured] we see, (a) and (d) (or (b) and (c)) differ only in the vertices labels, so they 
are topologically equivalent diagrams. These equivalent diagrams can be avoided easily by 
a simple trick. Let's look again with more care at the pairing procedure that leads Eq. (jSJ) 
to Eq. (jSJ). The external particle e\ is to be paired with a field if), and it faces two choices: 
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pairing with ipi in Hi or with ip 2 in H 2} that is 

{e 4 e 3 e 2 }T{[^ 1 (ei^i)</>i] [^^h]} 
{e 4 e 3 e 2 } T { [V'iV'i^i] [^ 2 (ei^ 2 )0 2 ] } 

Since ifi and H 2 are both "new" (i.e. not connected with any particles or fields), they are 
equivalent, thus the two choices in Eq. (J7J) are equivalent. Knowing this, we can suppress 
the equivalent diagrams by demanding e.\ to connect with Hi only, and by multiplying a 
factor 2 to take into account the other choice. The grouping for Eq. (jSJ) under this additional 
restriction becomes 

2 x IT { [MiXe^)] [(e 3 ^ 2 )(e 4 ^ 2 )(0i0 2 )] } 
2 x ^{^(e^)] [(e 2 V2)(e 4 V2)(0i0 2 )]} 

From Eq. (jSJ) we see, the equivalent diagrams disappear as expected. What's more, the 
factor 1/2! is canceled exactly by the multiplier 2. Indeed, this would happen for arbitrary 
1/N\ (only in the vacuum to vacuum case would this cancel be uncompleted, where a factor 
1/N remains.). 

The above example illustrates the first trick of our algorithm, namely 
Prescription 1: 

• In the procedure of an operator scanning over the fields in the T{} brace to pair with, 
when it encounter a "new" Hi, its scanning procedure is terminated after scanning 
over this Hi. The factor 1/N\ in Eq. ([3]) is discarded (or replaced by a 1/N in the 
vacuum to vacuum case) from the beginning. 

By this prescription, the equivalent Feynman diagrams due to vertex-relabeling are prohib- 
ited. 

Now, let's turn to another more complex example. We modify the previous example by 
change the physical model into 

H oc tfj e tf; e <f) n <f) n (9) 

This leads to 

S \2U = h {e4e 3 e 2 ei} x T { [0iV>i0i0i] [V> 2 ^ 2 2 2 ] } (10) 
Following the ordered pairing procedure and the prescription 1 above, we derive from Eq. 
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TO]) all the FD's as shown in Eq. (TTTj) and Figure [2j 



(11) 



T { [(e 2 $l)(e 1 V>l)(01a016)] [(e 3 ^2)(e 4 ^2)(02a026)] } 
r { [(e 2 ^l)(ei^l)] [(e 3 V ; 2)(e4V ; 2)(01a02a)(01b026)] } 
2 1 { [(e 2 V'l)(ei'0l)] [(e 3 ^ 2 ) (^2) (01602a) (01a02b)] } 
T { [(e 3 ^i)(ei?/'l)(01a01 O )] [(e 2 ^2)(e 4 V ; 2)(02a026)] } 
r { [(e 3 ^l)(eiV'l)] [(e 2 ^ 2 )(e 4 ^2)(01a02a)(01b02&)] } 
T { [(e 3 ^i)(ei^l)] [(e 2 ^2)(e 4 ^2) (01602a) (01a02b)] } 

We have add addition subscript a or b to the field in Eq. (jlip to distinguish different 




(d) 




FIG. 2: six Feynman diagrams corresponding to the six ways of combinations in Eq. (llli . 

pairings. Figure and Figure [2b are two unconnected diagrams, which should be rejected 
in real calculations, we leave this task to the programming stage. Here we care about the 
equivalent diagrams (b) and (c) (or (e) and (f)). The arising of these equivalent diagrams 
can be traced back to the pairing of 0i a with other fields. 0i a faces three choices: 0i o , 
02a and 02b, where the paring with 2a and with 2 (, are equivalent. Knowing this, we can 
constraint 0i a to pair with 2a only and represent the other equivalent choice by a multiplier 
2. With this additional restriction, and discarding the unconnected ones, the grouping for 
Eq. (10) becomes 



2 X T{[(e 2 ^i)(eiVl)] [(e 3 ^2)(e 4 ^ 2 )(01a02a)(01602o)] } 
2 X T{[(e 3 ^l)(eW>l)] [(e 2 ^2)(e 4 ^2)(01a02a)(01602 O )] } 



(12) 



This example shows the second trick of our algorithm, that is 
Prescription 2: 

• When an operator is paring with the fields in a vertex ifj, if there are m identical 
unpaired fields in Hi that can pair with this boson, the boson is restricted to pair with 
only one of them, and the other choices are represented by a multiplier m. 
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After the implementation of prescription 2, equivalent FD's originating from identical 
fields are greatly suppressed, and the symmetry factors are naturally generated. Note that 
the suppression by prescription 2 is incomplete, a small portion of equivalent FD's can still 
appear at two-loop order or higher. The correctness of the Feynman amplitude would not 
be hurt by the appearing of equivalent FD's, because the symmetry factor attached to each 
FD would naturally compensate for it, e.g. a (4*FD) may appear as two equivalent FD's of 
(1*FD) and (3*FD). We don't try to utilize a post eliminating procedure as done in many 
FD generators such as FeynArts, QGRAF and the grc part of GRACE to clean away the 
equivalent FD's, in order to keep the simplicity of the algorithm. 

In effective field theory one may encounter the case of multiple identical fermions appear 
in the interaction term, as 



In fact the two ^'s in Eq. (13) are not surely identical because of their different positions in 
H: when they are paired by the same fermion, a difference of a factor —1 might arise. For 
simplicity, in our algorithm we just treat these fermions as un-identical ones (see the input 
file description in section II III for more) . The readers might wish to find a more elegant way 
to deal with this problem. 

Minus signs coming from the interchange of fermions in the pairing procedure need also 
to be correctly produced, we leave this task to the programming stage. 

Overall, the "pairing order" the "prescription 1" and the "prescription 2" form the main 
parts of our algorithm for FD generation. 

III. A EXAMPLE C PROGRAM REALIZING THE ALGORITHM 

In this section, we briefly discuss a small C implementation of the algorithm proposed in 
section UT1 



The C code can be download freely at the website [16]. In the C code, two independent 
functions, judgeConnectionO and f ermionAcrossSignO , are made to perform the tasks 
of judging the connectivity of a FD and the sign from fermion interchanging respectively. 
judgeConnectionO is implemented after a FD is generated, f ermionAcrossSignO is 
implemented in the FD making stage: every time an fermion operator pairs with another 
fermion field, the f ermionAcrossSignO is called once. Other functions can be learned from 



(13) 
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the source code [ly] directly. 

The input file of the program is simple and clear. The following is an input file example 
corresponding to Eq. (10). 

The physical model (i.e. the relevent interactions) 

instructions : 

'U' means a Boson having no charge, i.e. it has no antiparticle , e.g. photon 
'C means a Boson which has antiparticle, e.g. W+ and W- 
'F' means a Fermion 

total number of interactions >>> 
1 

list of the interactions >>> 
3> Fl F-l Ul:2 



The physical process 



total number of incoming and outgoing particles >>> 
4 

the incoming and outgoing particles >>> 

2> Fl F-l 

2> Fl F-l 

total loops >>> 

1 

In the model-defining part of the input file, the £ Ui' or 'Ci' can be followed by a ' :n' which 
defines the multiplicity of the particle 'Ui' or 'Ci' in an interaction. If a user also want to 
write some interaction with multiple identical fermions, she should write as "...Fi...Fi..." , 
i.e. write each and every fermion explicitely. 

The outputs corresponding to the above input file look like: 

The physical model (i.e. the relevent interactions) 

>>> the interactions 
II: (F1,F-1,U1:2) 



The physical process 

>>> the external particles 



PI l 


(F-l) 


P2: 


(Fl) 


P3: 


(Fl) 


P4: 


(F-l) 


>» 


the i 


HI 




H2 




H3 




H4 





The Feynman diagrams 



(i) 

PI— >H1 (Il.leg[2]) 

P2— >H1 (Il.leg[l]) 

P3— >H2 (II. leg [1]) 

P4— >H2 (II. leg [2]) 

HI— >H2 (Il.leg[3]— >Il.leg[3]) 

HI— >H2 (Il.leg[3]— >Il.leg[3]) 

symmetryFactor = 2 * (1/2 ! ) * (1/2 ! ) = 1/2 

f ermionSign=-l 



(2) 

PI— >H1 (Il.leg[2]) 

P2— >H2 (II. leg [1]) 

P3— >H1 (II. leg [1]) 

P4— >H2 (Il.leg[2]) 

HI— >H2 (II. leg [3] — >I1. leg [3]) 

HI— >H2 (Il.leg[3]— >Il.leg[3]) 

symmetryFactor = 2 * (1/2 ! ) * (1/2 ! ) = 1/2 

f ermionSign=l 
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The C program is powerful, yet small. It is a totally general purpose FD generator, that 
is, it can receive arbitrary user defined models and arbitrary processes and any loops. Yet, 
benefit from the simplicity of the algorithm, the C program is small, i.e., ~ 500 lines of code 
as a whole, and ~ 200 lines as the core part for realizing the algorithm. It is also very fast: 
3860 FD's of uu — > ti at two loop order in QCD model can be generated in 0.015 seconds 
in a normal PC (compared to minutes by FeynArts). 

However, as only an illustration of our algorithm, it still lacks many functionalities which 
one might find useful: showing the Feynman diagrams in a graph view mode, picking out 
the one-particle irreducible diagrams, to name a few. 

IV. CHECKING THE CORRECTNESS OF THE ALGORITHM/C-PROGRAM 

In order to check the correctness of the proposed algorithm and the C-program, we use 
FeynArts and our C-program to generate FD's for the process uu — > ti in three different 
models 

MODEL1: uug + tig 
MODEL2: uug + tig + g 3 
MODEL3: uug + tig + g 4 

and at three different loop orders: tree, one-loop and two-loop. Table 1 shows the total 
number of FD's generated by the two programs, from which one see, the two programs 
agree with each other except for the MODEL2 and MODEL3 at two-loop order. These 
differences are attributed to the fact that our algorithm generates more equivalent FD's in 
these cases. In table 2 and table 3 we show explicitly all the equivalent FD's generated by 
our C-program and by FeynArts (When doing the comparison, we found that FeynArts also 
generate some equivalent FD's). Combining Table [H Table HT1 and table IHIl one see, the two 
programs agree exactly with each other. 

TABLE I: Comparison between the total numbers of FD's generated by our C-program and Fey- 
nArts. The numbers at the left (right) of the comma are from our C-program (FeynArts). 





MODEL1 


MODEL2 


MODEL3 


tree 
one-loop 
two-loop 


1, 1 
18, 18 
303, 303 


1, 1 
28, 28 
825, 764 


1, 1 
19, 19 
369, 355 
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TABLE II: Equivalent FD's generated by our C-program (first row) and FeynArts (second row) 
under the model MODEL2. The numbers in bracket are the ID number of the FD's generated in 
our C-program. 

(12)(15), (13)(16), (14)(17), (19)(20), (32)(35), (33)(36), (34) (37), (39)(40), (51)(53), 
(52)(60), (56)(61), (64) (67), (65)(66), (68)(70), (69) (72), (73) (74), (75)(78), (76) (77), 
(79)(81), (82)(97), (83)(98), (84)(99), (85)(100), (86)(102), (87)(103), (88)(104), (89)(105), 
(90)(107), (91)(108), (92)(109), (93)(110), (94)(111), (95)(112), (96)(113), (101)(121), 
(106)(122), (114)(117), (115)(118), (116)(119), (123)(126), (124)(125), (137)(139), 
(138)(146), (142)(147), (188)(190), (189)(197), (193)(198), (329)(332), (330)(333), 
(331)(334), (339) (340), (342) (345), (343) (346), (344) (347), (349)(350), (375) (376), 
(387) (389), (388) (396), (392) (397), (434) (436), (435) (443), (439) (444), (475) (478), 
(476) (479), (477) (480), (483) (484) in total: 66 pairs 

u I 




T177P2N731 T177P3N732 T178 P2 N738 T178P3N739 T179 P2 N745 T179P3N746 




T180P2N752 T180P3N753 T181P2N759 T181P3N760 intotaL 5pairS 



TABLE III: Same as Table U but with the model changed to MODEL3. 



(10)(11), (21)(22), (25)(26), (30)(31), (32)(33), (34)(35), (43) (45), (44) (46), (69)(71), 
(70)(72), (141)(142), (152)(154), (153)(155), (176)(178), (177)(179) in total: 15 pairs 



u t 

t tt 

H89P2N353 T189P3N354 in total: 1 pairs 



V. CONCLUSION AND DISCUSSION 

An algorithm for FD generation is proposed in this paper. It is simple in concept, easy for 
coding. A C-program realizing this algorithm is presented. It is small in size and run fast, 
yet is a totally general FD generator: it receives arbitrary user defined model and arbitrary 
process as input and generates FD's at any order. 

We hope this simple algorithm could make life easier for high energy researchers, who 
would like to make their own FD generators and add functionalities to their own taste. 
Researchers who are developing new calculation techniques may find the algorithm of special 
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interests, as it could help in both testing new ideas and constructing new calculation tools. 
It can also be useful for researchers working on the effective field theory due to its 'arbitrary 
user defined model' property. 
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